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ABSTRACT 

Magnetized neutron stars constitute a special class of compact objects harbouring 
gravitational fields that deviate strongly from the Newtonian weak field limit. More¬ 
over strong electromagnetic fields anchored into the star give rise to non-linear correc¬ 
tions to Maxwell equations described by quantum electrodynamics (QED). Electro¬ 
magnetic fields close to or above the critical value of Bq = 4.4 x 10® T are probably 
present in some pulsars and for most of the magnetars. To account properly for emis¬ 
sion emanating from the neutron star surface like for instance thermal radiation and 
its polarization properties, it is important to include general relativistic (GR) effects 
simultaneously with non-linear electrodynamics. This can be achieved through a 3-1-1 
formalism known in general relativity and that incorporates QED perturbations to 
Maxwell equations. Starting from the lowest order corrections to the Lagrangian for 
the electromagnetic field, as given for instance by Born-Infeld or Euler-Heisenberg 
theory, we derive the non-linear Maxwell equations in general relativity including 
quantum vacuum effects. We also derive a prescription for the force-free limit and 
show that these equations can be solved with classical finite volume methods for hy¬ 
perbolic conservation laws. It is therefore straightforward to include general relativity 
and quantum electrodynamics in the description of neutron star magnetospheres by 
using standard classical numerical techniques borrowed from Maxwell and Newton 
theory. As an application, we show that spin-down luminosity corrections associated 
to QED effects are negligible with respect to GR corrections. 

Key words: gravitation - magnetic fields - plasmas - stars: neutron - methods: 
analytical 


1 INTRODUCTION 

As a final stage of the stellar evolution process, neutron stars 
constitute a special class of compact objects showing strong 
field effects in both gravitational and electromagnetic inter¬ 
actions. Typical values for their radius lye around R~ 12 km 
whereas their fiducial mass is about M = 1.4 Mq (where 
Mq is the mass of the sun) thus leading to a compactness 
parameter defined by 

H = 2GM/Rc^ « 0.34 . (1) 

This compactness is close to the maximal reachable com¬ 
pactness defining a black hole. G is the gravitational con¬ 
stant and c the speed of light. General relativity is therefore 
required for an accurate description of the gravitational field 
in the vicinity of the neutron star. Moreover, from simple es¬ 
timates of their spin-down luminosity being interpreted as 
magneto-dipole losses, normal pulsars can harbour magnetic 
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fields of strength comparable lo B ~ 10® T and even higher 
fields for magnetars, of the order B « 10^°“^^^ T. About 
thirty magnetars ar e known today and com piled in a cat¬ 
alogue described in lOlausen fc KaspH (l2014l l . These values 
are comparable or largely above the critical field value of 

Bq = « 4.4 X 10® T. (2) 

me is the mass of an electron, e the absolute value of its elec¬ 
tric charge and ft the reduced Planck constant. At such field 
strengths, quantum electrodynamical corrections to classi¬ 
cal electromagnetism as given by Maxwell equations become 
significant. Vacuum polarization and electron-positron pair 
creation are two examples of such QED effects. The per¬ 
turbations induced by QED can be conveniently expressed 
in terms of an effective Lagrangian field theo ry such as the 
one found empirically bv iBorn fc Infeldl l|l934h or derived di¬ 
rectly from a first order expansi o n of t he QED Lagrangian 
as given by iHeisenberg fc Euleil (ll936l L Such effective La- 
grangians are very fruitful in getting more physical insight 
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and intuition on QED effects and represent very efficient 
tools to extend classical electrodynamics methods to the 
realm of quantum physics. Extreme fields around Bq can¬ 
not be reach in terrestrial laboratories although the power of 
current lasers are approaching it i n order to do reproducibl e 
experiments and check the theory (iKing fc Di Piaz^l2014h . 
Pulsars and magnetars could be used as extraterrestrial 
laboratories to investigate QED in strong electromagnetic 
fields and in curved space-time. The electromagnetic prop¬ 
erties of quantum vacuum in sp e cial relativity are su mma- 
rized in iBattesti fc Rizz3 (1201^ 1. iRufhni et al.l ll2010ll pro¬ 
pose a comprehensive review on electron-positron pair cre¬ 
ation/annihilation in the physical and astrophysical con¬ 
texts. The classical picture of high-energy radiation pro¬ 
cesses in astrophysics need also to b e revised according to 
QED corrections. This is discusses in Iffarding &: Tail (l2006l l 
with a special emphas ize to neutron star interiors and atmo - 
spheres. For instance. iHevl fc HernauistI (|2005 I 'I: Ih^ (120071 1 
suggested an explanation of the non-thermal emission from 
magnetar in terms of QED. 

Not only radiation is affected by QED but also the 
behaviour of pair plasmas and electr on-ion plasmas as ex¬ 
plaine d in depth in the review by lUzdenskv fc RightlevI 
ll2014l l. The plasma dielectric tensor and therefore also the 
normal modes of propagation of electromagnetic waves and 
plasma oscillations are affected. This had led some authors 
to extend the equations of fluid hydrodynamics t o some 
quantum aspects. For instance iManfredi Sz Ha^ (1200 il l 
looked at some of these corrections and iHaad 1 I 2 OO 5 I 1 added 
the magnetic field into the description leading to quantum 
plasma magnetohydrodynamic (QMH D). QED corrections 
to MH D have already been proposed bvlThompson fc Blaed 
dlOOSll and by Iflevl fc Hernauis^ (Il999l l . Interestingly, the 
non-linear nature of electrodynamics in strong fields induces 
shocks as in the case of flows in hydrodyna mics as pointed 
out b y Iflevl fc HernauistI lll998h . see also iMazur fc Hevll 
(I 2 OIIII for non-linear wave propagation in the magneto- 
spheric plasma. From a more fundamental point of view, 
QED corrections can be brought to any magnetic multipole, 
like for instance t he QED corrections to the m agnetic dipole 
as investigated bv lHevl fc HernauistI (Il997bll . 

Gravi t ationa l effects are usually ignored. Neverthe¬ 
less Ih^ (I 200 H 1 applied quantum field theory near a 
rotating black hole and showed tha t it could co p iously 
produce electron-positron pairs. Also iRufBni et al.l (I 2 OI 3 I 1 
studied QED corrections in strong gravity by inspecting 
spherically symmetric black holes u sing what they calle d 
the Einstein-Euler-Heisenb erg theory. IPenisov et al.l ll2004h ; 
iDenisov fc SvertilovI (l2005l l looked at electromagnetic wave 
propagation in a strong gravitational field described by the 
Schwarzschild m etric and including QED eff e cts. In another 
series of papers, IDenisov fc SvertilovI (l2003l ): IDenisov et al.l 
(|20141 examined the propagation of light rays, bending and 
time delay due to vacuum birefringence. 

From an observational point of view, X-ray polariza¬ 
tion is the key observable to study QED effects in gen¬ 
eral relativity in the context of neutron star magnetospheres 
and atmosphere s. Such signatur e s have indeed already been 
investigated by iHevl fc ShavivI (l2002f) . They showed that 
the X-ray polarization d egree is largely enhan ced compared 
to classical field theory. iHevl fc ShavivI (I 2 OOC 1 I ) also demon¬ 
strated a possible phase lag between different wavelengths. 


iTaverna et al.l (l2014l l showed that phase-resolved polarime- 
try gives insight into magnetar magnetospheres in the ultra 
strong field regime already with modest X-ray telescopes. 
At optical and infrared wavelengths, plasma and vacuum 
polarization effects compete, leading to possible constraints 
on the tota l charge density within t he magnetosphere as ex¬ 
plained bv IShannon fc Hevll (l2006ll . Moreover because vac¬ 
uum birefringence in strong magnetic fie ld leads to vari¬ 
able refractive indexes. IShaviv et al.l (Il999ll predicted strong 
QED lensing analogue to the gravitational lensing although 
in QED lensing the strength of th e lens depends on the po¬ 
lariza tion states of the photons (iDupavs et al.l I2OO5I : iKimI 
I2OI2II . Such effects are supposed to help diagnosing the mag¬ 
netic field structure around neutron stars. 

To study the magnetosphere of neutron stars, it is often 
useful to start with a simple approximation called force-free 
electrodynamics (FEE). It represents a zero order approx¬ 
imation to compute the plasma response and current den¬ 
sity knowing the electromagnetic field topology. FEE has 
been used in the last decade by many au thors to investi- 


gate pulsar and magnetar magnetospheres (iTimokhin 20061: 

Snitkovskv 

20061: McKinnev 20061: Petril 

2012 

: Parfrev et al.1 

2 OI 2 I: IPetri 

l2015al'). iFrevtsis fc Grallal ( 

2015 

) showed that 


quantum force-free electrodynamics can be recast into the 
3-1-1 language as a set of time evolution equations for the 
electric and magnetic field. The system looks very similar to 
special-relativistic FEE except for some corrections in the 
current density. 

In the birefringent medium induced by vacuum polar¬ 
ization and in the presence of plasma, vacuum resonances 
occurs, which translates i nto a mode convers ion from low to 
high opacity, according to iLai fc h 3 (l2003bh . But this reso¬ 
nance depends on the energy E of the photon, so it is im¬ 
portant to include possible gravitational redshift effects into 
that picture, especially when the radiation is emitted from 
the neutron star surface. This vacuum resonance imprints a 
special signature in the X-ray polari zation properties from 
its surface emission dhai fc Hdl2003al l. Because the vacuum 
resonance condition includes E it has to be corrected for the 
strong gravitational field. A 3-|-l formalism as the one we 
developed here might help to better understand such res¬ 
onance in neutron star magnetospheres. For instance, our 
framework could be used to extend the Mont e Carlo simula¬ 
tions for Compton scattering performed bv iBulik fc Mill^ 
(Il997h for soft gamma re peaters in flat space-time or by 
iFernandez fc Da^ (I2OIIII in a Schwarzschild background 
metric acc o untin g for light ray deflection as proposed by 
iHevl et ahl (l2003ll . It could also serve to compute the radia¬ 
tion spectra from vac uum polariz ation and proton cyclotron 
resonances as done bv lOzell (l2003ll or atmospheres of neutron 
stars and their resulting spectra as explained bv lHo fc Tail 
(I2OO3II . As a general conclusion, our approach combining 
gravity and QED in a 3-|-l split of space-time will be ex¬ 
tremely useful to straightforwardly extend the investigations 
cited before into the realm of general relativity. This will en¬ 
able to make accurate predictions about the observational 
signature of neutron star atmospheres and magnetospheres. 

In this paper, we include the most general background 
metric into the QED description. The 3-1-1 formalism used to 
split the covariant Maxwell equations into space and time is 
summarized in section (2] Next the Lagrangian formalism for 
the electromagnetic field is defined in curved space-time in 
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section [3l In section |4] we derive the non-linear electromag¬ 
netic field equations in curved space-time. An application 
to the spin-down luminosity in strongly magnetized neutron 
stars is discussed in section [5] Concluding remarks are given 
in section [6] 


2 THE 3+1 FORMALISM IN GENERAL 
RELATIVITY 


In order to apply traditional finite volume schemes and 
to get more physical insight into electromagnetism in gen¬ 
eral relativity, we remind for completeness the 3+1 lan¬ 
guage often used to transform covariant equations into time- 
dependent hyperbolic systems of spatial vectors in three di¬ 
mensions. To this aim, we split the four dimensional space- 
time into a 3+1 foliation such that the metric can be ex¬ 
pressed as 

ds^ = c? dF — Jab (da;“ + P°‘ cdt) {dx^ + cdt) (3) 


where x* = (ct,a;“), t is the time coordinate or universal 
time and some associated space coordinates. We use the 
Landau-Lifschitz convention for the metric signature given 
by (+,—,—,—) llLandau k, Lifchi^IlfiSfibl b a is the lapse 
function, /3“ the shift vector and jab the spatial metric of 
absolute space. By convention, latine letters from a to h 
are used for the components of vectors in absolute space 
(in the range {1,2,3}) whereas latine letters starting from 
i are used for four dimensional vectors and tensors (in the 
range {0, 1,2,3}). Our deriva tion of the 3+1 equation s fol¬ 
lows the method o utlined by Komlssaroy|j [2004l . l201ll 'l and 
extensively used bv iPetril Il20lll20l4l2015all . 

Let F'’*^ and be the electromagnetic tensor and 

its dual respectively. It is useful to introduce the following 
spatial vectors (B, E, D, H) such that 


(4a) 

Ea = ^eabaC*F^^ (4b) 

D°' = eocaF‘^° (4c) 

Ha^-^CabcF*’^ (4d) 

2/To 


where eo is the vacuum permittivity and /tq the vacuum 
permeability, Cabc = y/j £abc. the fully antisymmetric spatial 
tensor and Sabc the three dimensional Levi-Civita symbol. 
The contravariant analogue is j is the de¬ 

terminant of the spatial metric jab- The three dimensional 
vector fields are not independent, they are related by two 
important constitutive relations, namely 


£oE = QD + eoc/3 x B 

(5a) 

XT X. X D 

fiQ H = aH -. 

(5b) 


£o c 


The curvature of absolute space is taken into account by the 
lapse function factor a in the first term on the right-hand 
side and the frame dragging effect is included in the second 
term, the cross-product between the shift vector P and the 
fields. 

To complete the description of the electromagnetic field 
in general relativity including QED corrections, we need to 
derive the field equations. This is done starting from a La- 
grangian for the electromagnetic field, as presented in the 
next section. 


3 LAGRANGIAN OF THE 

ELECTROMAGNETIC FIELD 


In classical field theory, the Lagrangian C of the electromag¬ 
netic field is given according to the two field invariants ex¬ 
pressed in covariant form as X\ = Fik F''^ and X 2 = Fik 
To the lo west order in these invariants, the classical La¬ 
grangian l|jacksonll200ll l is given by 

Co = --^FikF^'^-FAi ( 6 ) 

4/^0 

where P is the four current and Ai the four potential of 
the electromagnetic field. In general relativity, the evolution 
equations for the electromagnetic field follow then from the 
variational principle exp ressed by Euler-Lagrange equations 
JUzan k. Deruellal2014r ) such that 


dc 1 ^ ^ dc 

dAi ddkAi 


( 7 ) 


where g = ce y/j is the determinant of the space-time metric 
in equation ®. Any non-linear covariant theory of electro¬ 
dynamics starts from a Lagrangian including higher orders 
of the invaria nts Xi and X 2 , like fo r insta nce t he Lagrangian 
propo sed by iHeisenberg fc Eul^ Jlfififil l or iBorn fc InfeTdl 
(Il934h . To remain as general as possible in our discussion, 
we use a parametrized post-Maxwellian description, simi¬ 
lar to the parametrized post-Newtonian case used for the 
gravitational field. The lowest order correction includes two 
parameters ( 771 , 772 ) such that the Lagrangian to this order 
becomes 


£. — To + 771 Xi + 772 X 2 . (8) 


In the Euler-Heisenberg prescription we have 


Ogf 1 

180 n 2 fio Bq 



(9a) 

(9b) 


whereas for the Born-Infeld Lagrangian in the weak field 
limit we have 


771 = 


1 

32 go 


(10a) 


772 = 771 (10b) 

with Psf the fine structure constant and h = 9.18 x 10^^ T 
the empirical maximal absolute field strength in Born-Infeld 
theory. 

Comparing the relative strength of the classical and 
QED Lagrangian, we deduce that quantum corrections to C 
are of the order 


C-Co 

To 


4/70 771 Fife 


Qsf 

90 7t Bq 


2.6x10” 


_ 

R2 


( 11 ) 


and therefore remain small even in the case of magnetic fields 
close to the critical field Bq . QED can always be treated as a 
small non-linear perturbation of Maxwell equation. For mag¬ 
netic fields well above Bq the aforementioned perturbative 
Lagrangians do not hold any longer but it can be shown that 
the relative strength between clas sical and QED Lag r angian 
remains much less than unity dHevl fc Hernauist Il997al : 
iRufiini fc Xuel[2006h unless B reaches unrealistic value of 
the order BQe^"'^“=f « 10®™ T because of the logarithmic 
dependence in ln(B/ BQ) of the Lagrangian in t his ultra high 
magnetic field limit llLandau fc Lifchit j|l989al . 
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4 NON LINEAR ELECTRODYNAMICS IN 
GENERAL RELATIVITY 

From the above first order corrections to the Lagrangian C, 
eq. m, it is straightforward to get the time evolution of the 
electromagnetic field, that is the non-linear Maxwell equa¬ 
tions in general relativity. 


4.1 Maxwell equations 

The equation of motion for the fields are given by the ex¬ 
pression 0. To write them down in the 3-1-1 language we 
note that 

o ^ Tplm 

+ +8 7]2*i^ X2 (12) 


/^O 

or briefly with the spatial vectors D and B 

dC _ F*"* 

ddlAm fJ-Q 

V £0 


^2 = 32 rj2 


D B 

£0 C 


(13a) 

(13b) 

(13c) 


Following the standard 3-1-1 decomposition used in general 
relativity we introduce two more auxiliary vector fields de¬ 
noted by E and G and defined by 


F ^ D + ^ B 

c 

G = (iH-^E . 


(14a) 


Straightforward computation shows that 

dC 

dAi 

■pOk 


= -r 


(15c) 


^ V F 

= (15d) 

V-g odkAa V go / 

= —!^ft(V7F)-lvxG. 
a ^7 a 

(15e) 

The inhomogeneous Maxwell equations then follow immedi¬ 
ately as 


It is seen from eq. GHl) that the primary fields to be evolved 
are B and F. The other auxiliary fields are then deduced 
from the constitutive relations. Indeed, B and F being 
known, we can retrieve D from eq. (I14all . Next from D and 
B we can get E and H through eq. 0. Finally G is ob¬ 
tained from eq. (I14bll knowing E and H from the previous 
calculation. This completes one full time step to advance the 
primary fields B and F. 

The source terms p and J are left free so far. In quantum 
vacuum as well as in classical vacuum, they vanish p = 0 
and J = 0. In the most general case, source terms have 
to be specified by some other equations, like the conserva¬ 
tion of energy-momentum of the plasma. Nevertheless, in 
some restricted problems, it is possible to compute the cur¬ 
rent density from only the knowledge of the electromagnetic 
field. Such an example, called force-free electrodynamics is 
presented in the next section. 


4.2 Force-free quantum electrodynamics 
(FFQED) 

The source terms have not yet been specified. Having in 
mind to apply the above equations to pulsar and magnetar 
magnetospheres, we give expressions for the current density 
in the limit of a force-free plasma, neglecting inertia and 
pressure. The force-free condition in covariant form and in¬ 
cluding QED corrections reads 


(14b) 

Fife /'= = 0 

(17) 

and in the 3-1-1 formalism it becomes 



J ■ E = 0 

(18a) 

(15a) 

pE4-JxB = 0 

(18b) 


which implies E • B = 0 and therefore also D • B = 0. Thus 
the parameter ^2 must vanish and therefore it follows that 
E ■ B = 0. Applying the usual technique from special rela¬ 
tivistic electrodynamics we get the current density including 
general relativity and QED such that 


^ E A B BVxG — FVxE^ 
J = p - B. 


52 


R2 


(19) 


Because = *F“''nfc and D“/eo = F“''nfc, B and D/eo 
can be interpreted as the magnetic and electric field respec¬ 
tively as measured by the fiducial observer whose four ve¬ 
locity is rife = (a c, 0) . Moreover 


V • F = p (16a) 

V X G ^ J +-^dti^F) . (16b) 

V7 

Note that p = a P/c has to be interpreted as the external 
charge density and J“ — a F as the external current den¬ 
sity. Vacuum polarization is treated implicitly through the 
definition of the vectors F and G as given by equations d. 
The homogeneous Maxwell equations are not modified, they 
read as in classical general relativity 


V 


xE = —^dt{^B) 
v/7 


V ■ B = 0 . 


(16c) 

(16d) 


Pnk=pP ( 20 ) 

thus p is the electric charge density as measured by this 
same observer. Its electric current density j is given by 

aj = J + PC/3. (21) 

Maxwell equations m, the constitutive relations 0,111 
and the prescription for the source terms set the background 
system to be solved for any prescribed metric in the low elec¬ 
tromagnetic field limit according to the classical Lagrangian 
correction. As a check of these equations, let us compute 
their simplified version in the classical general-relativistic as 
well as in the special-relativistic QED limits. 
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4.3 Classical general-relativistic limit 

If QED corrections are neglected, we have to set ^1 = 1 
and ^2 = 0. In this limit the new auxiliary fields F and 
G reduce respectively to the classical fields D and H. The 
inhomogeneous Maxwell equations become 

V D = p (22a) 

VxH = J+-^dt{fi7D) (22b) 

V7 

which are the expre ssions found bv iKomissarovI (l2004l . [201lh 
and bv IPetril (l2013l h The current density then simplifies to 

, ExB BVxH — DVxE„ , , 
J = +- ^2 -B (23) 

as expected from the above cited works. 


4.4 Quantum special-relativistic limit 


In the other limit, when the gravitational field becomes neg¬ 
ligible, the lapse function is equal to one, a = 1, and the shift 
vector vanishes, (3 = 0. The constitutive relations ([5ll sim¬ 
plify to So E = D and po H = B. Maxwell equations then 
resemble those in a medium given by 


V X E = 

V xG = 

with the initial conditions on 
B fields such as 


dB 

dt 


J 



the divergence of the 


(24a) 
(24b) 
F and 


V • B = 0 (24c) 

V ■ F = p. (24d) 


To retrieve the more familiar notation, we should substi¬ 
tute F hy D and G by H. Electrodynamics in the pres¬ 
ence of strong electromagnetic fields can be described by 
non linear Maxwell equations derived from an effective La- 
grangian computed by Euler and Heisenberg. Quantum elec¬ 
trodynamics described vacuum as a polarized and magne¬ 
tized media without external current density, J = 0, and 
without any charge density, p = 0. The quantum vacuum is 
depicted by a magnetization M and a polarisation P such 
that 


F = So E + P 

(25a) 

G = B/po-M . 

(25b) 

To the first order in the fine structure constant, the 
and Heisenberg Lagrangian shows that 

Euler 

P = k,{2[E‘^ - c B'^) E + 7 c {E ■ B) B) 

(26a) 

M = -k{2 F (E^ -F B^)B-7 F {E ■ B) E) 

(26b) 

with 


Oisf 

45 7t pQ d'- Bq 

(27) 

Consequently, our constitutive relations include both limits, 
the classical general-relativistic field and QED in Newtonian 
gravitational held. We refer to it as general-relativistic force- 
free quantum electrodynamics (GRFFQED). 


5 SPIN-DOWN LUMINOSITY OF PULSARS 
AND MAGNETARS 


What kind of change in the spin-down luminosity can we ex¬ 
pect from the vacu um polarization? We know a lready from 
previo us works by iRezzolla fc AhmedovI (l2004h and IPetril 
(l2013ll that general relativity leads to an increase by a fac¬ 
tor 2 to 6 of the inferred spin-down losses compared to flat 
space-time. Nevertheless Maxwell equations remain linear 
in general relativity. The enhanced luminosity can be inter¬ 
preted as a combination of magnetic field amplification and 
gravitational redshift of the angular frequency of the neu¬ 
tron star due to the curvature of space-time. Moreover, from 
the lowest order corrections brought by QED, Maxwell field 
equations become non-linear. These non-linearities will per¬ 
turb the topology of the dipole but they will also generate 
higher order multipoles like for instance an hexapole due to 
terms containing products of (E^, B^) with {B, E). For the 
corrections to the dipole we have 


(JEdip ^ 4asf ^ ^ 

Edip “ B ~ 457t B2 ~ ^ B2 


(28) 


which remain s negligible. These conclu sions have already 
been found bv lHevl fc HernauistI (Il997bh . For the hexapole, 
using the point multipol e formula for the most luminous 
mode i.e. m = 2 given bv IPetril (l2015bli we find 


Ehex ^ ^ (Rfi\ ^ ^ ^ — 

Edip 25 Vruy* B^ ^ B^ ^ B^ 


(29) 


so even less significant than the dipole corrections. We 
conclude that whatever the period and the magnetic field 
strength of pulsars and magnetars, QED corrections to the 
spin-down luminosities of any multipolar component remain 
meaningless with respect to the corrections brought by gen¬ 
eral relativity alone. 


6 CONCLUSION 

In this paper we showed how to include vacuum polarization 
into the description of a force-free magnetosphere in general 
relativity according to the weak electromagnetic field limit 
i.e. in the first order perturbation theory of the Lagrangian 
for the electromagnetic field. By introducing two new aux¬ 
iliary fields, it is possible to cast the full set of Maxwell 
equations into a classical three dimensional picture treating 
curved space-time and vacuum polarization as two medium 
with specified constitutive relations for both parts. Our new 
formalism should help to quantify the merit of both contri¬ 
butions, general relativity and quantum electrodynamics, to 
the dynamics of neutron star magnetospheres, especially for 
magnetars and pulsars with high-B fields. As shown in the 
previous section, on a global length scale related to the size 
of the magnetosphere and to its rotational braking, QED 
is irrelevant to account for variation in the spin-down lu¬ 
minosity. Although QED corrections to the electrodynamics 
of neutron star magnetosphere remain weak, its implications 
for interpretation of observations like propagation and polar¬ 
ization of electromagnetic waves and pair creation are likely 
to be important. 
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